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Abstract — This paper presents theoretical and experimental 
results showing the effects of linear and circular birefringence 
on the sensing coil of an in-line Sagnac interferometer current 
sensor. Linear birefringence induces local scale factor variations 
along the length of the sensing coil which can be suppressed with 
circular birefringence. We show that by winding the sensing coil 
with tension in a helix configuration, these variations may be 
effectively reduced to zero for ordinary single mode fiber. In our 
opinion, this method of implementation overcomes the need to 
use annealed fiber in the sensing coil. 

Index Terms — Fiber-optic current sensors, Sagnac interferom- 
eter. 



I. Introduction 

FIBER-OPTIC current sensors which rely on Ampere's 
law and the Faraday effect are universally sensitive to 
changes in the linear birefringence, 8, in their sensing coils. 
Changing stresses and temperature may change the total 6 or 
its distribution along the sensing coil which in turn changes 
the evolution of, the state of polarization of the light through 
the sensing coil. This in turn changes the point by point 
sensitivity of the coil to the magnetic field. The result is that 
both the scale factor relating the sensor output to the current 
in the enclosed wire and the isolation of the sensor to currents 
in nearby wires are subject to errors and changes in real 
environments [l] f [2 J. In this paper, we address these errors 
particularly as they affect the in-line Sagnac interferometer 
current sensor [3]-[8] shown in Fig. 1. We present theoretical 
and corroborating experimental results showing how the point 
by point sensitivity of the sensing fiber varies along its length 
in the case of a uniformly distributed 8. By integrating over 
the local sensitivity function, we obtain the overall sensor scale 
factor. Results for both open and closed-loop operation of the 
sensor are obtained. 

A particular concern is to find a method to desensitize the 
system against the effects of 6. The use of annealed fibers 
[9] or specialty low-strain optic coefficient fibers [10] solves 
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Fig. I. In-line Sagnac interferometer current sensor. 

the problem, however, it is desirable in practice to avoid 
these brittle fibers both because of the added handling cost 
and the loss of the significant advantage the in-line Sagnac 
interferometer version of the current sensor possesses which 
is that it can be installed without disconnecting the current 
carrying wire by simply wrapping the sensing coil around 
the conductor. We present theoretical and experimental results 
that explicitly show how adding circular birefringence to the 
sensing coil stabilizes the scale factor on a point by point 
basis. Winding the sensing coil under modest tension in a loop 
helix is demonstrated to add sufficient circular birefringence to 
an ordinary single-mode fiber to provide for a constant scale 
factor along the entire length of the sensing fiber. 

The scale factor of the in-line Sagnac interferometer is 
determined both by the quality of the quarter-waveplate and 
the sensing fiber. For mitigating the errors associated with the 
use of an imperfect quarter waveplate, we refer the reader to 
our previous paper (8]. 

n. In-Line Sagnac Interferometer Current Sensor 

The in-line Sagnac interferometer current sensor (Fig. I) is 
a reciprocal device which extracts the Faraday induced phase 
shift by interfering two orthogonally polarized light waves 
traveling through the optical circuit in reverse order. These 
waves travel down the polarization maintaining (PM) fiber 
delay line linearly polarized and through the sensing head 
circularly polarized. On the return trip, each of the optical 
waves ideally assumes the orthogonal sense of polarization 
at each point in the sensing head and PM fiber delay line. 
Reciprocal optical effects arising from strains induced by 
acoustic vibrations and changing thermal environments are 
cancelled because both beams travel the same optical path 
[3], [6]. For current sensing applications, the only contributor 
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to the- phase difference between the two beams is that from 
the Faraday effect, and the total phase shift is 4VNI where 
V is the Verdet constant of the sensing fiber glass, N is the 
number of turns of fiber around the current carrying wire, and 
/ is the current in the wire. For simplicity, we substitute F 
for ViVJ, i.e., F == VNI. 

m. Mathematical Modeling of the 
Local Scale Factor 

We analyze the spatial variation of the scale factor relating 
the response of the sensor to a localized magnetic field (Fig.. 2). 
We assume that the total linear and circular birefringences, (6 
and 2T), are uniformly distributed throughout the sensing coil 
of length L. 

In an optical fiber current sensor, the sensing coil is nor- 
mally wound around the current carrying conductor to utilize 
Ampere's law in extracting the information of the current. 
In this case, the magnetic field H is integrated all along the 
loop which consists of an integral number of turns around the 
conductor. Here, we are interested in the variation of scale 
factor along the path of integration, namely, along the sensing 
coil. To model this, we use a Jones matrix analysis assuming a 
localized magnetic field fi to be located a distance z from the 
terminating mirror of the sensing coil. The following parameter 
definitions are used to simplify the notation: 



-y=(J)v^ 2 + (27T 

P = «.(!) 

* = sin(f) 
*=cos(l) 

y = sin (|) cos x 

Z = sin(2) 



* cos x 



•sin* 



siux- 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 

(7) 
(8) 
(9) 



The light first travels through the polarizer, L pt after which it 
is aligned to a 45° splice, £45°. to the polarization maintaining 
(PM) fiber delay line LpMm which possesses a delay time of 
r/2 and a differential polarization delay time of At. A bire- 
fringence modulator located just before the PM fiber delay line 
is included within this Jones matrix. This modulator provides 



the phase difference modulation between the two interfering 
light waves. As in the fiber optic gyro, the modulator operates 
at maximum efficiency when the round-trip delay time, r, 
and the fundamental modulator frequency, / p , are related by 
f p = r/2. We modulate the birefringence using a waveform 
ip(t) provided by a piezoelectric tranducer (PZT — open-loop 
case) or an integrated optics modulator (IOM — closed-loop 
case). After propagating through the delay line, the light then 
passes through a quarter- waveplate, L X / 4 , aligned at 45° with 
respect to the delay line, £45°, before entering the sensing coil. 

The sensing coil is divided into three sections; two sections 
having lengths L — z and z which contain only the linear 
and circular birefringences surround a very short piece of 
fiber which is subjected to a magnetic field. The matrix 
L co \\\ represents the segment of fiber between the quarter- 
waveplate (A/4) and the concentrated magnetic field element, 
represented by L^in. The fiber element between the position 
of the magnetic field element and the mirror is represented by 
£coii2- We assume all elements are ideal with the exception 
of the sensing coil elements (L C oiji and L co ii2) which contain 
both 6 and 2T. 

The lightwave passes through L col ii and then through 
the concentrated magnetic field element containing Faraday 
rotation, Lr\n- It then continues through L C oii2. strikes the 
terminating mirror, L m i rro r. and retraces its way back through 
the optical circuit. To describe the propagation of the oppo- 
sitely traveling lightwave, we take the transpose of the each 
Jones matrix, reverse the sign of the off-diagonal elements, 
and further reverse the sign of the Faraday angle F. This 
convention preserves a right handed coordinate system with 
light propagation always to be in the +z direction. This 
convention accounts for the difference of the two Jones 
matrices depicting Faraday rotation in the sensing coil, Lrin 
and Z/Fout- The composite Jones matrix describing the total 
light propagation path is summarized in (20) 



-*PMout 
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^out - L p • L 4 $c • Lp MoiIt • L 45 o * L X /4 • L coin 

' ^Fout * ^««iI2 ' ^mirror ' ^roi!2 * ^Fin " ^coill 
* L\/ A * L450 " £pMin * £40* • L p • £? in . (20) 

From this output electric field, we calculate the total output 
light intensity striking the photo-detector of the sensor. Here, 
we assume a broadband light source and integrate over its 
frequency spectrum, 1/, thus eliminating all terms containing 
a cos2nisAr or sin2*-^r dependency. We also define the 
phase difference modulation, V m (t) t as 

* m («)MW^-r), (2!) 
The output intensity is given by 

/o»t = (Zw E* ut ) (22) 

Utilizing sinusoidal modulation, $ m (t) = ^ m cos^, to- 
gether with open-loop demodulation signal processing com- 
monly used for fiber-optic gyroscope applications, the first 
harmonic component, V lH9 of the sensor's output intensity 
signal is proportional to the Faraday rotation caused by a 
current 

^i/f = 2/ in . J,(* m ) 

* {[(4QY(-(PX) + RZ)(X 2 - V' 2 + Z 2 ))Jsin2F 
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Y 2 + Z 



2 ) 2 J sin4Fj 



(23) 

where.J in = (I^„| 2 >. 

We: assume the usual case that the Faraday induced phase 
shift, is small. sin2F % 2F and sin4F ^ 4F. We extract 
the birefringent dependent coefficient of the local scale factor 
by normalizing V lH to the light source intensity (/ in ), the 
modulation depth (J x (* m )), and the Faraday rotation (4F) 

Replacing the parameters P, Q. R, X, Y, and Z with their 
representations depicting 6 and 2T, (24) reduces to 

SF I<K . al (*, 2T. ,)=( m* + P«»(1:JS* + (2r)*)\ 
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Fig. 3. Scale factor along imperfect sensing coil, T - 0°, = 20°, and 
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(25) 



An inspection of (25) reveals that a ^independent local 
scale factor is a necessary and sufficient condition to ensure 
that linear birefringence is not affecting the total scale factor. 
For the case that the sensing coil contains no circular birefrin- 
gence (2T = 0), the local scale factor along the sensing fiber, 
SF| (M . al (£, varies as the top of a cosine curve (26) 

SF|oral(*, 2T = 0, z) = cos(r) ^ cos 6. (26) 

The local scale factor, including both linear and circular 
birefringence (S and 2T) in the sensing coil can be stabilized 
either by reducing S by annealing the sensing fiber or by 



0 0. 1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 I 
Fig. 5. Scale factor along imperfect sensing coil, T = 240°. fi = 20°, and 

adding a large amount of circular birefringence so as to swamp 
the linear birefringence. SF, ocaI (^2r, z) tends to unity for all 
z y provided 2T > 6. The fiber may be twisted to achieve 
a circular birefringence large enough to swamp the linear 
birefringence, but twisted fiber tends to untwist over time. 

It is interesting to note that as 2T increases, not only 
does the maximum deviation from unity decrease, but the 
periodicity of the variation also increases. Figs. 3-6 show the 
trend of scale factor variation along the sensing coil for varying 
amounts of IT for the two cases of an evenly distributed S 
of 20 and 40°. 

TV. Mathematical Modeling 
of the Overall Scale Factor 

A. Open-Loop Signal Processing Configuration 

Tn the open-loop version of the sensor, the birefringence 
modulator is driven with a sine wave modulation as described 
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above. The first harmonic portion of the signal at the pho- 
todetector is taken to be the system output. When the total 
Faraday induced phase shift is small, the total sensor response 
may be calculated by integrating the local response of the 
sensor over its length. Assuming as before that there exists 
uniformly distributed linear and circular birefringences, (25) 
may be integrated to obtain the birefringence dependent part 
of the overall scale factor SF(£, 2T). 



SF(*. 2T) = ± SF locaI (S, 2T, z) dz 



■( 



(2Tf + 6 2 cos(^6* + (2T) 2 ) 
6 2 + (2T) 2 

6 2 + (2T) 2 



(27) 



A particularly interesting open-loop signal processing 
scheme for the in-line Sagnac interferometer consists of 
dividing the first harmonic output by the peak intensity 
in the detected waveform. This scheme was shown in [8] 
to desensitize the sensor to imperfections in the quarter- 
waveplate. The peak output intensity of the current sensor 
scales with delta and 2T as 

W(*.2T) oc 1 + I S2 + (2T)2 I * 

(28) 

Performing the peak-divide signal processing scheme, the 
birefringent dependent coefficient of the Faraday rotation 
scales as dividing (27) by (28) 



SF P eak(^2T)oc 



SF(g,2T) 

Jpeak(£, 2T) 



(29) 



When using the peak divide signal and open-loop pro- 
cessing scheme, the total birefringent dependent scale fac- 
tor SF,,«>ak(^ 2T) is attenuated in the presence of linear 
birefringence, 6 t in the sensing coil. However, as circular 
birefringence, 2T is added, the scale factor error may be 
eliminated. Fig. 7 show the birefringent dependent scale factor 
as a function of 2T for 8 = 20° and 40°. As very large amounts 
of 2T are added, the scale factor error is effectively reduced 
to zero. 
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In the absence of circular birefringence in the sensing coil 
(2T = 0), the scale factor dependence, as a function of S t is 
given as 

SF(«, 2T = 0) = - J SF IocaI (£, 2T = 0, z) dz 



sin 28 
28 * 



(30) 



The dependence of the scale factor on linear birefringence 
when employing the peak divide signal processing scheme, 
SF P eak(<5, 2T = 0), is given below. Even in the absence 
of any linear birefringence compensation, the peak divide 
signal processing scheme does provide some scale factor error 
reduction (compare Figs. 7 and 8). 



SF pea k(<5 1 2T = 0)oc 



oc 



SF(g,gr = o) 

W(«,2T = 0) ~ l -fcos 2 * 



(31) 



Linear birefringence within the sensing coil near the quarter- 
waveplate (A/4) may simply be taken as an error in the 
A/4 itself, and is therefore compensated by the peak di- 
vide signal processing. However, linear birefringence existing 
further within the sensing coil yields an error that is not 
completely compensated by this signal processing. Thus the 
desirability of adding circular birefringence or eliminating the 
linear birefringence in the sensing coil remains. 

B. Closed-Loop Signal Processing Configuration 

Similar to the fiber-optic gyroscope, the in-line Sagnac 
interferometer current sensor may be operated "closed loop." 
A Ti-indiffused LiNbOa integrated optics straight waveguide 
modulator is used to provide the birefringence modulation. A 
phase ramp waveform is added to the usual periodic waveform 
for providing a non reciprocal phase shift equal and opposite 
to the Faraday induced phase shift. The first harmonic, Y*i// of 
the periodic waveform is thereby nulled in the sensor output. 

The Jones matrix modeling of the closed-loop sensor is the 
same as the open-loop sensor with the one exception that (12) 
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is changed to 

£PM out = e >(^ 



0 



(32) 



Here, 70 (5, 2T, F, z) represents the effective nonreciprocal 
phase shift that is added to the system by the birefringence 
modulator. After performing the Jones matrix analysis, the 
phase shift 70 which is required to zero V lH is given below 
as a function of linear birefringence, 6 t circular birefringence, 
27\ Faraday rotation F and position along the sensing coii z. 



7o(£,2T, F,z) = -2 tan" 1 { tan 2F 



■(■ 



/ (2T)* + * 2 cos(f + (2T)2) 
\ (2r) 2 + ^cos(v/6 2 + (2T) 2 ) 



))■ 



It is evident from (33) that a ^-independent local scale factor 
is again a necessary and sufficient condition to ensure that 
there is no adverse effect of linear birefringence on the sensor. 
As with the open-loop case, this condition is satisfied when 
2T » 6. 

To calculate the total phase shift which must be induced 
to zero V\n we perform a similar mathematical analysis 
which does not assume a localized magnetic field H % but 
rather assumes an evenly distributed magnetic field along 
the entire length of the sensing coil. For the closed-loop 
signal processing configuration, it is not appropriate simply to 
integrate the localized phase shift 7o(£, 2T, F t z) with respect 
to 2 to find the total phase shift 7 0 (£»27 , I F) required to 
zero V\h> The principle which allows the total scale factor 
to be derived by integrating the localized scale factor in the 
open-loop signal processing configuration is the fact that a 
linear sensor response is assumed. For the closed-loop case, 
which can handle much larger values of Faraday rotation, the 
assumption of a small Faraday angular rotation should not be 
pressed. 

The resultant total phase shift required to zero V lH in the 
closed-loop signal processing approach is given in (34), and 
the closed-loop birefringent dependent coefficient SF,-i(£,2T) 
given in (35) is shown in Fig. 8 „ 



7o(*,2T, F) = (-4F). 



(2T) 2 + 6 2 cos( vA 2 + (2T) 2 ) 



/ 
(34) 

(35) 



S W ,2T ) = 2^^. ' 

An interesting aspect of the closed-loop system is that the 
error has several zero crossings (Fig. 8) irrespective of the 
amount of total S in the sensing coil (assuming 2T > S). 
Precise and exact amounts of 2T may be introduced to any 
sensing coil, regardless of its given $ and the scale factor 
error associated with S can be completely overcome. The unity 
crossings of 70(^1 2T, F) (Fig. 8) are given below. 
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tan 2T = 2T. 



(36) 



However, as can be seen by Fig. 8, as the total amount of 
2T continues to increase, the error in the scale factor becomes 
negligible. A more robust solution is to induce as much circular 
birefringence as possible. 

V. Experimental Results 

A. Straight Sensing Coil 

We first performed an experimental analysis of the setup 
shown in Fig. 2. To apply a localized magnetic field H along 
the sensing coil, we used a very small air core inductor and 
mapped the magnitude of the scale factor as we moved the 
inductor along the sensing coil- For this series of tests, we 
simply twisted the fiber to induce 2T. The experimental results 
(Figs. 9-11) are seen to follow qualitatively the predicted 
behavior as shown in Figs. 3-6. 

Fig. 9 shows the scale factor response along the length of the 
sensing coil (one meter long) in the specific case of no circular 
birefringence, i.e., 2T = 0. Clearly the scale factor drops off 
as the top of a cosine curve. Fig. 10 shows the scale factor 
response along the length of the sensing coil when the entire 
length of the sensing coil has undergone a moderate amount of 
twisting. Its behavior resembles that shown in Figs. 4-5 for the 
oscillatory effect as well as decreased deviation from unity of 
the scale factor. Fig. 1 1 shows the scale factor response when 
a very large amount of twisting had been applied to the sensing 
coil. For this case, a large amount of 2T has been introduced. 
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and the scale factor is approximately constant all along the 
sensing coil. 

B. Loop Helix Sensing Coil 

An alternative way to add a large amount of circular 
birefringence to the sensing coil is to wind it in a helix [1 1], as 
shown in Fig. 12. This method of adding circular birefringence 
is much more robust than twisting the fiber, as the circular 
birefringence is the result of the geometrical torsion, and 
therefore does not change with temperature and lime. For our 
laboratory experiment testing this particular implementation 
of the sensing coil, our setup included a loop helix sensing 
coil wound with a sensing coil radius of approximately 30 
cm. Our helix contained twelve complete rotations evenly 
distributed around the sensing coil. The radius of our helix 
is approximately 5 mm. 

In this configuration, the sensing coil undergoes an incre- 
mental change of coordinate axes all along its length, provided 
there exist well defined axes of linear birefringence oriented 
parallel and perpendicular to the plane of the fiber bending. To 
ensure that the sensing coil contains well-defined fi axes which 
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follow the geometrical twist along the helix, the sensing coil 
should be held with some axial tension. Unless the sensing 
coil is very small, the bend birefringence is not typically large 
enough to overcome the intrinsic birefringence built into the 
fiber. However, bending with a small amount of axial strain 
does provide sufficient linear birefringence to define a rotating 
coordinate system. 

Figs. 13 and 14 show experimental confirmation tha.t wind- 
ing the sensing fiber in a helix and adding axial strain stabilizes 
the scale factor along the sensing loop. A localized magnetic 
field was induced using a small air-core inductor around the 
sensing fiber. When the helically wound fiber was loose, the 
cosine dependence of the scale factor along the length of the 
sensing fiber is evident as shown in Fig. 13. Then, when a 
small amount of tension was added to the helically wound 
fiber, the scale factor straightened out as shown in Fig. 14. 

VI. Conclusion 

We have performed a mathematical analysis of scale factor 
error associated with linear birefringence 6 in the sensing coil 
of the in-line Sagnac interferometer current sensor. In doing 
this, we found the spatial dependence of the scale factor along 
the sensing coil of the current sensor as well as the total scale 
factor of the current sensor as a function of both linear and 
circular birefringence (6 and 2T) in the sensing coil. 

We achieved experimental verification of Jhis spatial scale 
factor dependence along a sensing coil containing varying 
amounts of 5 and 2T. In doing so, we showed that the 
scale factor error of a sensing coil containing sufficiently 
large amounts of 2T may be effectively eliminated. We 
achieved a flat scale factor result experimentally, and we are 
therefore assured that our sensing coil is free from any linear 
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birefringence induced uncertainties. The use of a helically 
wound sensing coil under modest tension was shown to 
desensitize the sensor to intrinsic linear birefringence built into 
the sensing fiber. We believe this design obviates the need for 
the use of annealed fiber in the sensing coil. 
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